Pfaffian-like ground state for 3-body-hard-core bosons in ID lattices 
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We propose a Pfaffian-like Ansatz for the ground state of bosons subject to 3-body infinite re- 
pulsive interactions in a ID lattice. Our Ansatz consists of the symmetrization over all possible 
ways of distributing the particles in two identical Tonks-Girardeau gases. We support the quality of 
our Ansatz with numerical calculations and propose an experimental scheme based on mixtures of 
bosonic atoms and molecules in ID optical lattices in which this Pfaffian-like state could be realized. 
Our findings may open the way for the creation of non-abelian anyons in ID systems. 

PACS numbers: 03.75.Fi, 03.67.-a, 42.50.-p, 73.43.-f 
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Beyond bosons and fermions, and even in contrast 
to the fascinating abelian anyons (AA) non-abelian 
anyons (NAA) [2J exhibit an exotic statistical behavior: 
If two different exchanges are performed consecutively 
among identical NAA, the final state of the system will 
depend on the order in which the two exchanges were 
made. NAA appeared first in the context of the frac- 
tional quantum Hall effect (FQHE) , as elementary ex- 
citations of exotic states like the Pfaffian state 0, , the 
exact ground state of quantum Hall Hamiltonians with 
3-body contact interactions. Recently, the possibility of 
a fault tolerant quantum computation based on NAA 
has boosted the investigation of new models containing 
NAA |fj, as well as the search for techniques for their 
detection and manipulation JJ. Meanwhile, the versatile 
and highly controllable atomic gases in optical lattices 
Q have opened a door to the near future implementa- 
tion of those models as well as for the artificial creation 
of non- Abelian gauge potentials [!j . 

All actual models containing NAA are 2D models. The 
motivation of the present work is the foreseen possibility 
of creating NAA in one-dimension (ID). This long-term 
goal requires in the first place to define the concept of 
NAA, which is essentially 2D, in ID. For abelian anyons 
(AA) this generalization has been already made by Hal- 
danc Within his generalized definition the spinon 

excitations of ID Heisenberg antiferromagnets are classi- 
fied as |- AA ■ This classification becomes very natu- 
ral through the connection between the ID antiferromag- 
netic ground state (for a long-range interaction model, 
the Haldane-Shastry model ^H) and the Laughlin state 
|12| for bosons at v — 1/2. In a similar way we anticipate 
that a connection can be established between quantum 
Hall models containing NAA and certain long-range ID 
spin models exhibiting NAA within a generalized defini- 
tion 0. 

Here, far from analyzing the above questions in gen- 
eral, our aim is to pave the way for the creation of ex- 
otic Pfaffian-like states in ID systems, which we believe 
may serve as the basis to create NAA. We present a re- 
alistic ID system whose ground state is very close to 
a Pfaffian-like state. The actual system we consider is 
that of bosonic atoms in a ID lattice with infinite re- 



pulsive 3-body on-site interactions, which we call 3-hard- 
core bosons. Inspired by the form of the fractional quan- 
tum Hall Pfaffian state for bosons 0, ^| , we propose an 
Ansatz for the ground state of our system. This Ansatz 
is a symmetrization over all possible ways of distribut- 
ing the particles in two identical Tonks-Girardeau (T-G) 
gases [la , Il6) . Comparison of the Ansatz with numeri- 
cal calculations for lattices up to 40 sites yields very good 
agreement. As for fractional quantum Hall systems, NAA 
may be created here by creating pairs of quasiholes, each 
quasihole being in a different cluster Q . This possibility 
will be discussed elsewhere [lflj . 

Three-body collisions among single atoms rarely occur 
in nature. However, they can be effectively simulated by 
mixtures of bosonic particles and molecules. This has 
been proposed by Cooper |17| for a rapidly rotating gas 
of bosonic atoms and molecules. Here, we show that a 
system of atoms and molecules in a ID lattice can in 
a similar way effectively model 3-hard-core bosons. We 
will show that the conditions to realize this situation lie 
within current experimental possibilities. 

3-hard-core bosons. We consider a system of bosonic 
atoms in a ID lattice with repulsive 3-body on-site inter- 
actions. This system is described by the Hamiltonian: 



H 



(1) 



where the operator a\ (cti) creates (annihilates) a boson 
on site £, t is the tunneling probability amplitude, and 
U3 is the on-site interaction energy. From now on we 
will consider the limit U3 — > 00. In this limit the Hilbert 
space is projected onto the subspace of states with oc- 
cupation numbers ni = 0, 1,2 per site. We will refer to 
bosons subject to this condition as 3-hard-core bosons. 
The projected Hamiltonian has the form 



H* = 



(2) 



where the 3-hard-core bosonic operators 0,3^ obey 
( a 3,e) 3 — an d satisfy the commutation relations 
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These opera- 
tors can be represented by 3 x 3 matrices of the form 
/0 1 \ 

a 3 i = I V2 • In contrast to the usual hard- 
ly / 

core bosonic operators which are directly equiva- 
lent to spin- 1/2 operators, the operators Og, a 3 are re- 
lated to spin-1 operators {S + , S~ , S z } in a non-linear 

way: a 3 — > S + + — 1)S Z \ This mapping leads 

to a complicated equivalent spin Hamiltonian (with third 
and fourth order terms) which seems hard to solve. 

In the following we present an Ansatz wave function 
for the ground state of Hamiltonian (J2J). Our Ansatz is 
inspired by the form of the ground state for fractional 
quantum Hall bosons subject to three body interactions 
SO- [0|- The reason to believe that this inspiration may 
be good is the deep connection already demonstrated for 
the case of two-body interactions between ground states 
of certain ID models and those of 2D particles in the 
lowest Landau level (LLL) [Tl|. 

Let us now turn for a moment to the problem of bosons 
in the LLL subject to the 3-body interaction potential 



5 2 {zi ~ Zj)5 2 (zi - Zk) |3j, with z, = Xi+iyi being 



the complex coordinate in the 2D plane. For infinite in- 
teraction strength the exact ground state of the problem 
is the Pfaffian state [1[T^|: 



$ 3 cx S 




(3) 



This state is constructed in the following way. Particles 
are first arranged into two identical v = 1/2 Laughlin 

states .12|, $f cx Y\f<J( z t ~ z jf labeled by a =|,L 
Then the operator symmetrizes over the two "vir- 
tual" subsets of coordinates {zj} and {z\}. Note that the 
Laughlin state <E>2 °f each cluster is a zero energy eigen- 
state of the 2-body interaction potential Y^i^j $( z i ~ z i)- 
This guarantees that in a state of the form © three parti- 
cles can never coincide in the same position: for any trio, 
two of them will belong to the same group and cause the 
wave function $3 to vanish. 

In direct analogy with equation (0 wc propose the 
following Ansatz for the ground state of Hamiltonian J2Jl : 
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This Ansatz has the same structure as Q, but 
the Laughlin state has been substituted by a Tonks- 
Girardeau (T-G) state HI, cx Yiflj |sin« - arj)j, 
with x% = 2tt/Mi, i = 1,...,M, M being the num- 
ber of lattice sites. This state is the ground state of 
hard-core ID lattice bosons with Hamiltonian #2,0- = 

ditions fij 



£+1 + h.c.) and periodic boundary con- 
Here, a 2 fj are hard-core bosonic operators 



satisfying (a 2(T ) 2 = 0. Written in second quantization 

the Ansatz |@J takes the form: |* 3 ) = V ® l^zf), 

where V is a local operator of the form V = Vf M , and 
Ve is an operator mapping the single-site 4-dimensional 
Hilbert space of two species of hard-core bosons to the 
3-dimensional one of 3-hard-core bosons (see Fig. 1). 
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FIG. 1: Schematic representation of the operator Vi map- 
ping the single-site 4-dimensional Hilbert space of two species 
of hard-core bosons to the 3-dimensional one of 3-hard-core 
bosons. 

Let us analyze different characteristic properties of our 
Ansatz. Taking into account the well known result for a 
T-G gas, namely the scaling of the one-particle corre- 
lation function (a\ +A ae) as A -1 / 2 for large A 0], we 
can derive the following asymptotic behavior for the one- 
body and two-body correlation functions for the Ansatz 
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The proof of J5j is more involved and we will give just 
numerical evidence from our calculations bellow (see 
inset of Fig. 3). The two-body correlation jnj) is indeed 
in our case the (one-particle) correlation function for 
on-site pairs. This means that whereas the system seen 
as a whole exhibits some kind of coherence (the spatial 
correlation decaying slowly as A^ 1 / 4 ) the underlying 
system of on-site pairs is in a much more disordered state 
(with a fast correlation decay as A" 1 ). This is in contrast 
to what happens in a weakly interacting bosonic gas in 
which coherence between sites is independent of their 
occupation number. We can also obtain analytical ex- 
pressions for the relative occupation of single and doubly 
occupied sites. The average number of doubly occupied 
sites is ri2 = (alalaiai) /2 = z^ 2 /2, and the one of single 
occupied sites is given by n\ = {a\ai{2 — n^)) = v(2 — v). 
This distribution is very different from the Poissonian 
one, for which we have n\^° /n\° = vj2. 
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FIG. 2: Overlap (^exl^a). The main plot shows the depen- 
dence of the overlap with the filling factor v — N/M for a 
system of M = 10 (orange circles) and M — 20 (brown cir- 
cles) lattice sites. The inset shows the decrease of the overlap 
with increasing system size M, at v = 1. 



merically, the overlap is still good (« 0.90). 

Fig. 3 shows the statistical distribution of doubly and 
single occupied sites for l^ex), which is very close to the 
one of the Ansatz \$fa), and clearly different from the 
Poissonian distribution typical of a weakly interacting 
Bose gas. Fig. 4 shows the momentum distribution for 
particles and on-site pairs together with the long-range 
scaling of their spatial correlation functions. We can 
clearly see how the exact state exhibits all characteristic 
behaviors that we have discussed above for the Ansatz. 

Experimental proposal. 

Inspired by Cooper's ideas 01 f° r 2D rotating Bose 
gases we present an experimental scheme for the real- 
ization of Hamiltonian J2J. Let us consider a system of 
bosonic atoms and diatomic Feshbach molecules trapped 
in a ID optical lattice. The Hamiltonian of the system 
is H = H K + Hp + Hi Hill, where 



H K = -t a ^(aja i+ i + h.c.) - t m ^(mjm i+ i + h.c), 



As an additional property, the Ansatz has particle- 
hole symmetry. This means that for a filling factor of 
the form v = N/M = 2 — 77, with N the number of 
particles, the state we propose is just the Ansatz for holes 
at i>h — rj. However, as we can clearly see from its matrix 
representation, the Hamiltonian does not exhibit this 
symmetry. This tells us that our Ansatz may not work, 
as we will see, for the whole regime of filling factors. 

Numerical calculations. In order to determine the qual- 
ity of our Ansatz we have performed a numerical cal- 
culation for the ground state of Hamiltonian To 
obtain the numerical ground state |^ e x) we have used 
variational Matrix Product States (MPS) of the form 

£?, SA ,Tr(,i; : ....IV; -,...*v; 0, with matri- 

ces A of dimension D = 12(15), and d = 3. We can 
estimate the error of this calculation to be smaller than 
10~ 5 for the system sizes (M < 40) we have considered. 
To calculate the overlap of |^ex) with our Ansatz ^3) 
we first construct the MPS state that best approximates 
|\l/3) for a given D. This is done in the following way. We 
first build the MPS ground state of Hamiltonian H^. We 
then take the tensor product of this state with itself, ob- 



taining a MPS with d = 4, which is closest to 



The dimension of the matrices of this state is very large 
and we use a reduction algorithm [2(| to reduce it to the 
initial size. Finally we apply the operator V by local ten- 
sor contraction and normalize the resulting d = 3 MPS 
state. For the matrix dimensions we used the error made 
was always smaller than 10~ 3 . 

The results are shown in Fig. 2. The main plot shows 
the overlap (vE'cxI 1 ^) as a function of the filling factor 
v = N/M for a fixed system size. We find very good 
overlaps (0.98-0.96) for v < 1.25. For v > 1.25 the over- 
lap decreases. The inset shows the overlap as a function 
of increasing system size M, at fixed filling factor v = 1. 
At M — 38, the maximum size we have considered nu- 



H F = 



H, 



EA t , Uaa ft 1 9 , f 1 l \ 

Amimj + —^-a\a\aiai + —^=(771^04 + h.c), 



A I 



m\a\a % mi + ^ m\m\m,imi . (7) 



A I 



i = l 



Here, the bosonic operators for atoms (molecules) a; (mi) 
obey the usual canonical commutation relations. The 
Hamiltonian Hk describes the tunneling processes of 
atoms and molecules, occurring with amplitude t a and 
t m , respectively. The term Hp is the Feshbach reso- 
nance term, with A being the energy off-set between 
open and closed channels, U aa the on-site atom-atom 
interaction and g the coupling strength to the closed 
channel. Hamiltonian Hi describes the on-site atom- 
molecule and molecule-molecule interactions. We will 
assume a situation in which U aa ,U a m,U mm > 0, and 
A > 0. Furthermore, we will consider the limit in which 
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FIG. 3: Average occupation of sites with one (left figure) and 
two particles (right figure) for the exact (solid lines), Ansatz 
(dotted lines) and Poissonian distribution (dot-dashed lines), 
as a function of the filling factor v. The system size is M — 20. 
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FIG. 4: Quasi-momentum distribution of particles oc 
E(A e ~' tA ( a f+A fl ') (orange, left axis) and of on-site pairs 
oc J2 e A e~' kA (al +A al +A agae) (blue, right axis), with 
k — 2n/Mn, n = 0, . . . , M — 1. Results are shown both for the 
exact ground state (solid lines) and the Ansatz (dashed lines). 
The inset shows the long-distance scaling of the correlation 
functions {a\ +A a t ) ~ A - " 1 (orange), and {al +A a\ +A a e a e ) ~ 
A"" 2 (blue), for the exact ground state (circles, ct\ = 0.22, 
a.2 = 0.83) and the Ansatz (diamonds, a\ = 0.24, o?2 = 0.99. 
Parameters are M = 20 and v = 1. 

7 2 = g 2 /2A 2 <C 1. Within this limit the formation of 
molecules is highly suppressed due to the high energy 
offset, A. However, virtual processes in which two atoms 
on the same lattice site go to the bound state, form a 
molecule and separate again, give rise to an effective 3- 
body interacting atomic Hamiltonian of the form p9|: 

i i 

- ^ 7 2 ^(( a l) 2 K +1 ) 2 +h.c.) 

i 

+ (U aa -g 2 /A)J2(4) 2 M 2 , (8) 

i 

where we have neglected higher order terms in 7 2 . As- 



suming U aa — g 2 /A, and t m -f 2 <C t a , H c ff reduces to 
Hamiltonian Q with t = t a and C/3 = U am l 2 - Finally, 
assuming U am ^f 2 ^ t a we end up with the desired Hamil- 
tonian @ for 3-hard-core bosons. 

Let us now summarize the requirements and approxi- 
mations we have imposed and discuss their experimental 
feasibility in typical setups with 87 Rb. W e have assumed 
,9 2 /A = U aa . Since g = v / U aa AnAB/2 21], with AB 
being the width of the Feshbach resonance and A/i the 
difference in magnetic momenta, we need A = AfiAB /2. 
For the Feshbach resonance at 1007. AG, this implies 
A/h = 441kHz [23L l24l|. Furthermore, we have assumed 
7 2 < 1, t m -i 2 < t a and U am ~f 2 /t a > 1. Written in 
terms of the lattice and atomic and molecule param- 
eters we have 7 2 = (a 3 a ^a/a 2 ± )(T]E R / A/jAB), 
Ua,n/t a = (V6/4tt) {a^a / o 2 )^ 2 exp(+7r 2 ?? 2 /4) and 
t m /t a = 2exp(-7rV/4), where X] = {Vq/Er) 1 / 4 , with 
Vo the lattice depth, Er = h 2 /(8ma 2 ) the recoil en- 
ergy, m the atomic mass, and a the lattice constant. 
The parameters a 3 ® and 0^ are the 3D scattering 
length fo r atom-a tom and atom-molecule collisions, and 
a± = y/ h/uo±_ra is the transversal confinement width, 
with lo± the transversal trapping frequency. Assuming 
typical values rf = 50(70), a™ ~ 5nm |lf, a?£ ~ a 3 ° 
|2r| . a = 425nm, and = 27r x 20kHz [23, we ob- 
tain: 7 2 = 3.55 (3.86) x 10" 3 , t m /t a = 5.29 (0.22) x 10~ 8 , 
and U am "/ 2 /ta — 1.35(30.4) x 10 3 , clearly satisfying the 
required conditions. 

Regarding detection of the Pfaffian-like ground state, 
the characteristic difference between both the momentum 
distribution and number statistics of particles and on-site 
pairs could be observed via spin-changing collisions 

In conclusion, we have shown that the ground state of 
3-hard-core bosons in a ID lattice can be well described 
by a Pfaffian-like state which is a cluster of two T-G 
gases. We have shown that such a state may be accessi- 
ble with current technology with atoms and molecules in 
optical lattices. We believe that our findings may open a 
new path for the creation of NAA. 

B. Paredes would like to thank M. Greiter for estimu- 
lating discussions and critical reading of this manuscript. 
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